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ABSTRACT 


A  method  is  developed  for  calculating  the  growth  of  a  turbulent 
boundary  layer  at  hypersonic  Mach  numbers.  Excellent  agreement  with 
experimental  results  from  axisymmetric  nozzles  has  been  obtained  by 
the  application  of  this  method.  The  method  utilizes  a  modification  of 
Stewartson's  transformation  to  simplify  the  integration  of  the  momentum 
equation.  Heat  transfer  is  taken  into  account  by  evaluating  the  gas  prop¬ 
erties  at  Eckert's  reference  temperature  and  by  using  a  modification 
of  Crocco’s  quadratic  for  the  temperature  distribution  in  the  boundary 
layer.  A  new  empirical  relation  is  used  for  the  incompressible  friction 
coefficient  which  agrees  with  experimental  data  over  a  Reynolds  num¬ 
ber  range  from  10®  to  10’. 


NOMENCLATURE 


A  Flow  area 

a  Speed  of  sound 

Cf  Local  skin-friction  coefficient 

Cf  Mean  skin-friction  coefficient 

h  Enthalpy 

H  Boundary  layer  form  factor 

M  Mach  number 

N  Parameter  in  power-law  equations  for  skin  friction 

P  Function  defined  in  Eq.  (4) 

P  Pressure 

0  Function  defined  in  Eq.  (4) 

R*  Radius  of  curvature  of  nozzle  at  the  throat 
Reynolds  number  based  on  x 
Rx  Incompressible  Reynolds  number  based  on  X 
R(9  Reynolds  number  based  on  6 

r  Local  distance  from  axis  of  flow  to  wall 

s  Distance  measured  along  wall  surface 

T  Temperature 

U  Transformed  velocity  component  parallel  to  wall 
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u  Velocity  component  parallel  to  wall 

X  Transformed  coordinate  along  axis  of  flow 

X  Coordinate  along  axis  of  flow 

Y  Transformed  coordinate  normal  to  wall 

y  Coordinate  normal  to  wall 

a.  Function  defined  in  Eq.  (20) 

|3  Function  defined  in  Eq.  (21) 

y  Ratio  of  specific  heats 

5*  Boundary-layer  displacement  thickness 

A  Boundary-layer  thickness 

6  Boundary-layer  momentum  thickness 

H  Dynamic  viscosity 

P  Density 

Tyi  Shearing  stress  at  wall 

CO  Angle  of  wall  surface  with  respect  to  axis 

Subscripts 

aw  Value  at  adiabatic  wall  conditions 

e  Value  at  free- stream  static  conditions 

i  Value  for  incompressible  flow 

o  Value  at  free-stream  stagnation  conditions 

s  Local  stagnation  temperature 

t  Total 

tr  Transformed 

w  Evaluated  at  the  wall 

Superscripts 

*  Evaluated  at  M  =  1  except  5* 

Evaluated  at  reference  temperature 
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INTRODUCTION 


Many  investigators  have  studied  the  problem  of  calculating  the 
growth  of  turbulent  boundary  layers.  For  the  incompressible  adiabatic 
case,  such  calculations  are  fairly  straightforward  since  they  depend 
only  upon  the  empirical  relationship  of  the  skin-friction  coefficient  with 
Reynolds  number  when  the  pressure  gradients  in  the  direction  of  flow 
are  favorable  or  absent  (Ref.  1).  Even  in  the  presence  of  adverse 
pressure  gradients,  empirical  methods  have  been  devised  for  the 
purpose  of  predicting  the  location  of  separation  (Ref.  2). 

As  the  speed  of  the  gas  outside  the  boundary  layer  increases,  the 
effects  of  compressibility  on  the  skin-friction  coefficient  must  also  be 
tEiken  into  account.  For  many  cases  in  supersonic  flow,  the  heat 
transfer  between  the  gas  and  the  wall  may  be  neglected  and  methods 
such  as  Tucker's  (Ref.  3)  may  be  used  to  calculate  the  boundary-layer 
growth.  There  is  enough  difference,  however,  between  the  free- 
stream  static  temperature  and  the  adiabatic  wall  temperature  that  some 
intermediate  "reference"  temperature  must  be  used  for  evaluating  the 
compressible  skin-friction  coefficient  in  order  to  obtain  good 
correlation  with  experimental  results.  Tucker  uses  a  reference 
temperature  equal  to  the  arithmetic  average  of  the  adiabatic  wall 
temperature  and  the  free-stream  static  temperature. 

For  gases  with  a  Prandtl  number  of  1,  the  adiabatic  wall  tempera¬ 
ture  is  equal  to  the  stagnation  temperature  and  the  thermal  boundary 
layer  has  the  same  thickness  as  the  velocity  boundary  layer.  For 
gases  such  as  air  which  have  a  Prandtl  number  less  than  1,  the 
adiabatic  wall  temperature  is  less  than  the  stagnation  temperature 
and  the  thermal  boundary  layer  is  thicker  than  the  velocity  boundary 
layer.  Bartz  (Ref.  4)  attempts  to  take  this  effect  into  account, 
together  with  the  effect  of  heat  transfer,  by  obtaining  simultaneous 
solutions  of  the  momentum  equation  and  the  energy  equation.  He, 
however,  restricts  his  solution  to  a  one -seventh-power  velocity  profile 
and  uses  a  reference  temperature  equal  to  the  arithmetic  average  of 
the  wall  temperature  and  free-stream  static  temperature. 

When  the  wall  temperature  is  much  lower  than  the  adiabatic  wall 
temperature  (often  the  situation  at  hypersonic  speeds),  the  reference 
temperature  should  be  somewhat  higher  than  the  arithmetic  mean 
between  the  wall  and  free-stream  static  temperature  as  shown  by 
Sommer  and  Short  (Ref.  5)  and  Tendeland  (Ref.  6).  The  reference 
temperature  suggested  by  Eckert  (Ref.  7)  correlates  with  existing 
data  as  well  as  that  used  by  Sommer  and  Short.  Persh  (Ref.  8) 
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attempts  to  account  for  this  same  effect  by  using  a  reference  tempera¬ 
ture  equal  to  the  temperature  at  the  edge  of  the  laminar  sub -layer, 
where  this  temperature  is  defined  as  a  function  of  velocity  ratio 
according  to  Crocco's  quadratic  modified  by  using  the  adiabatic  wall 
temperature  to  account  for  the  Prandtl  number  being  less  than  one. 

One  important  use  of  an  accurate  method  of  calculating  turbulent¬ 
boundary-layer  growth  at  hypersonic  Mach  numbers  is  in  the  design  of 
axisymmetric  hypersonic  wind  tunnels.  The  perfect -fluid  (potential- 
flow)  contour  cam  be  accurately  calculated  by  the  method  of  character¬ 
istics,  but  the  attainment  of  uniform  flow  in  the  test  region  depends 
also  upon  the  accuracy  of  the  boundary-layer  correction  to  the 
theoretical  contour.  Each  of  the  available  methods  of  calculating 
boundary-layer  growth  possessed  some  drawback.  This  fact  led  to  the 
development  of  the  method  described  herein,  which  utilizes  a  modifi¬ 
cation  of  Stewartson's  transformation  (Ref.  9)  to  aid  in  the  integration 
of  the  momentum  equation.  The  main  difference,  however,  between 
this  method  and  other  methods  is  in  the  evaluation  of  the  compressible 
skin-friction  coefficient  and  the  transformed  form  factor.  In  order  to 
simplify  the  calculations,  a  new  empirical  equation  was  developed  for 
the  incompressible  skin-friction  coefficient  which  agrees  with  experi¬ 
mental  data  over  a  range  of  Reynolds  number  (Rx;)  from  10'  to  10*. 
Excellent  correlations  have  been  obtained  between  experimental  values 
of  boundary-layer  displacement  thickness  and  those  calculated  by  this 
method. 

Since  this  method  was  developed  primarily  for  hypersonic  wind 
tunnels,  it  may  not  be  suitable  where  adverse  pressure  gradients  are 
present  or  where  the  temperatures  are  sufficiently  high  to  cause 
dissociation. 


DEVELOPMENT  OF  METHOD 


MOMENTUM  EQUATION 


For  the  purpose  of  this  report,  steady  compressible  flow  is 
assumed.  Thus,  the  von  Karman  momentum  equation  for  axisymmetric 
flow  can  be  written  in  the  form: 


d0 

ds 


2  -  M* 

+  H 

dM  1 

dr 

f  y  ~ 

1  ,  \ 

ds  r 

ds 

—  M’  j 

- 

(1) 


In  order  to  obtain  a  solution  for  the  momentum  thickness  the  values 
of  M,  11,  r,  ajid  Cf  must  be  known  as  a  function  of  distance  s  along  the 
surface.  For  many  cases,  such  as  wind  tunnel  nozzles,  these  values 
are  known  as  a  function  of  distance  x  along  the  center  line.  Since 
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ds 

dx 


sec  CO 


(2) 


Eq.  (1)  may  be  written  as 


d^  ^  ^ 

+ 

1 

CM 

1 _ ^ 

dM  1  1 

dr"'  _ 

dx 

+  m^) 

dx  r 

dx 

2 

sec  CO 


When  CO  is  small,  the  assumption  that  sec  co  =  1  can  be  made  with 
negligible  consequences. 


(3) 


It  can  be  recognized  that  Eq.  (3)  is  a  linear,  first-order,  ordinary 
differential  equation  of  the  form 


+  d  P(x)  =  Q(x) 
dx 


(4) 


which  has  the  solution 

0  (x)  =  e  l^/Q  (x)  e  dx  +  constant  J 

Since  P(x)  and  Q(x)  are  generally  non-analytic  functions  of  x,  numerical 
methods  must  be  used  to  evaluate  the  indicated  integrations. 

The  solution  of  Eq.  (3)  is  considerably  simplified,  however, 
through  the  use  of  the  equations 

0  (l  +  =  0tr  (6) 


and 


H  =  Htr(l  +  M^)  + 


(7a) 


or 


H  +  1  =  (Htr  +  1)  ^  (71- 

which  are  obtained  from  Stewartson's  transformation  as  shown  in 
Appendix  A.  By  substituting  Eq.  (6)  and  (7)  into  Eq.  (3),  the  equation 


dx 


^  (2  + 
dx 


H») 


^tr  dr 
r  dx 


Cf 

— i-  sec 
2 


y+  1 


(8) 
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is  obtained.  By  multiplying  Eq.  (8)  by  rM  ^  ,  the  left-hand  side 

becomes  a  perfect  differential  if  Htr  is  assumed  to  be  constant  over 
the  interval  of  integration.  Therefore, 


M 


y+i 

£l  2(y-i) 

2  V  r  Q 


dx 


(9) 


Thus,  the  solution  of  the  momentum  equation  is  reduced  to  a  single 
integration  which,  in  general,  must  be  accomplished  by  numerical 
methods. 


Substitution  of  Eq.  (6)  into  Eq.  (9)  yields 


y+  1 

y  +  1  ■ 

0rM 

b 

2  +  Htr 


sec 


w 


y  +  1 

2  “  \  / 


dx 


(10) 


It  may  be  noted  that,  as  long  as  Hu  as  defined  by  Eq.  (7)  is  approxi¬ 
mately  constant,  Eq.  (6)  and  (7)  may  be  considered  to  be  merely 
integrating  adds  irrespective  of  the  validity  of  Stewartson’s  transfor¬ 
mation. 


When  there  is  no  heat  transfer  and  the  Prandtl  number  is  one,  the 
transformed  form  factor  becomes  equal  to  the  incompressible  value. 

In  an  attempt  to  account  for  heat  transfer  and  for  Prandtl  numbers  other 
than  one,  use  is  made  of  Crocco's  quadratic  for  the  temperature 
distribution  in  the  boundary  layer  according  to  the  equation 


X 

Te 


(11) 


which  is  the  form  used  by  Persh  (Ref.  8).  As  shown  in  Appendix  A,  the 
use  of  this  equation  yields 


Combining  Eqs,  (7)  and  (12)  gives 


(12) 


Hu 


(13) 


which  relates  the  transformed  form  factor  in  terms  of  the  incompressible 
value  and  the  temperature  ratios. 
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For  the  evaluation  of  Eq.  (1)  for  a  particular  problem,  values  of 
r,  M,  cd,  and  y  are  given  as  functions  of  x  .  With  the  wall  temperature 
given  or  assumed  as  a  function  of  x,  Htr  may  be  found  by  Eq.  (13)  in 
terms  of  Hi.  Still  remaining  to  be  determined  are  Cf  and  II j,  which  are 
considered  in  subsequent  sections. 


DETERMINATION  OF  Cf 


The  compressible  skin-friction  coefficient  Cf  in  Eq.  (10)  is  a 
function  of  Mach  number,  Reynolds  number,  and  heat  transfer  and  is 
defined  in  terms  of  the  shearing  stress  at  the  wall. 


2  Pe  (14) 

In  the  evaluation  of  Cf,  efforts  are  made  to  obtain  empirical  correlations 
with  the  incompressible  value  which  is  a  function  of  Reynolds  number 
based  either  on  the  momentum  thickness  or  distance. 


F(H0.)  =  G(Rxi) 


(15) 


The  assumption  has  been  made  that  there  is  some  reference  point  in 
the  compressible  boundary  layer  where  the  temperature  and  density 
are  such  that 


where 


fji: 

2 


=  Fine')  =  G(ax') 

p 


(16) 


R0'  =  ^ 

P 

(17) 

Rx'  = 

(18) 

P 

and  the  values  of  p'  and  p'  are  evaluated  at  the  reference  temperature. 
Therefore, 


C,  _  p'  Cf  ^  Je  C(' 

2  Pe  2  T'  2 


(19) 


since  the  pressure  is  assumed  to  be  constant  through  the  boundary 
layer. 


It  may  be  shown,  however,  that 

F(R0')  ^  yV  G(Rx') 
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by  using  the  power-law  equations  for  skin-friction  coefficient 

F  (R0')  =  - 

(Ri9')  ^ 

and 


(20) 


where 


G  (Rx')  = 


A 


1 _ 

h  N  +  1 


/3 


N  +  1 


1 

/n±_L\  n  +  1 

\  N  J 


(21) 


(22) 


can  be  derived  from  the  incompressible  relationships.  In  the  absence  of 
pressure  gradient,  Eq.  (1)  can  be  written  as 


<i6  _  A. 

dx  2 


li 


(23) 


then 


2 


A  f(R0')  = 


(24) 


(25) 


and,  after  integrating, 


e 

X 


JS__ 

N  +  1 


N  +  1 


On  the  other  hand,  if 


G(Rx') 


T' 


_ ft _ 

I 

\  IL  ' 


(26) 


(27) 
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then 


de 


_ _ 

X 


and,  after  integrating, 


T 

^  e 

t' 


R 

_ N _ 


_J _ 

N  +  1 


Te 

T' 


f_N  +. 

L.  _N 


a 


\ 

N  +  1 


1 

N  +  1 


(28) 


(29) 


It  is  obvious  that  Eq,  (26)  is  not  equal  to  Eq.  (29),  and  therefore 

Eq.  (24)  and  Eq.  (27)  cannot  both  be  correct.  Experimental  correlation 

must  be  made  to  determine  which  equation  gives  the  better  results. 

The  determination  of  the  proper  value  of  reference  temperature  to 
be  used  must  also  depend  upon  experimental  correlation.  Eckert 
(Ref.  7)  suggests  the  relation 


T'  =  0.5  Tw  +  0.22  Taw  +  0.28  Te 
which,  if  the  recovery  factor  is  equal  to  0.  896,  becomes 


(30) 


T'  =  0.5  Tw  +  (0.5  +  0.0394  M')Te 
Sommer  and  Short  (Ref.  5)  suggest  the  relation 

T'  =  0.45  Tw  +  (0.55  +  0.035  M')  Te 

which,  again  if  the  recovery  factor  is  equal  to  0.  896,  becomes 


(31) 


(32) 


T'  =  0.45  Tw  +  0.195  Taw  +  0.355  Te 


(33) 


Eckert's  relation  shows  slightly  greater  effects  of  both  heat  transfer 
and  Mach  number  than  does  that  of  Sommer  and  Short  and  appears  to 
give  slightly  better  correlation  with  the  meager  amount  of  data  which 
exists  at  hypersonic  Mach  numbers.  In  order  to  take  into  account  the 
effects  of  variable  specific  heat,  Eckert  suggests  the  use  of  a  reference 
enthalpy 


h  —  0.5  hw  0.22  h^w  0,28  he 


(34) 


for  evaluating  the  physical  properties  of  the  gas.  Although  the 
derivations  described  herein  are  developed  for  constant  specific  heat, 
this  modification  can  be  easily  substituted  throughout. 
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The  correlation  of  theoretical  results  with  experimental  data  can 
be  made  by  comparing  the  ratios  of  compressible  skin-friction 
coefficients  to  the  incompressible  value  at  the  same  free-stream 
Reynolds  number.  Most  of  the  experimental  data  have  been  obtained 
in  wind  tunnels.  In  the  supersonic  range  up  to  a  Mach  number  of  about 
5,  stagnation  temperatures  are  of  the  order  of  100°F  or  560°  R  while 
the  static  temperatures  decrease  with  increasing  Mach  number  until  a 
value  of  about  100°  R  is  reached.  At  higher  Mach  numbers,  the 
stagnation  temperature  is  increased  to  maintain  the  static  temperature 
at  about  100°R  in  order  to  avoid  liquefaction  of  the  constituents  of  the 
air.  The  temperature  variation  which  is  suggested  for  correlating 
wind  tunnel  results  is  shown  in  Fig.  1  along  with  Eckert's  reference 
temperature  for  the  adiabatic  wall  condition.  The  temperature  range 
is  so  great  that  a  simple  power  law  for  the  variation  of  viscosity  with 
temperature  is  not  valid,  and  Sutherland's  law  must  be  used  above 
198.7°R.  Below  this  temperature,  there  is  some  meager  evidence 
that  the  viscosity  varies  proportionally  with  the  temperature  (the  curve 
is  tangent  to  the  Sutherland  curve).  Misleading  results  can  be  obtained 
if  such  factors  are  not  taken  into  account  for  correlation  over  a  wide 
range  in  Mach  numbers. 

A  correlation  of  Cf/Cf.  as  a  function  of  Mach  number  is  shown  in 
Fig.  2.  The  experimental  data  were  taken  from  Refs.  10  through  17 
and  are  for  approximately  adiabatic  wall  conditions.  The  values  of  Cf. 
used  for  the  ratios  were  obtained  from  the  relation 

0.088  (log  -  2. .3686) 

'i  =  OogHi,^  1.5)^  (35) 

which  is  developed  in  Appendix  B.  The  values  of  Cf  were  obtained 
from  the  relation 


Te  0.088  (log  Rx  '  -  2-3 686) 

t'  (log  R,,'  -  1.5)'  (36) 


where 


P 


TV' 


(37) 


and  llx,  the  free-stream  Reynolds  number,  is  the  same  as  R*;  for  this 
curve.  The  static  and  reference  temperatures  from  Fig.  1  were  used 
to  obtain  the  ratio  of  Cf/Cf.  shown  in  Fig.  2.  The  good  correlation 
shown  indicates  that  Eq.  (36)  provides  an  adequate  method  of  determining 
Cf  for  use  in  Eq.  (10).  It  should  be  noted  that  Eq.  (36)  is  of  the  general 
form  of  Eq.  (27)  which  therefore  represents  experimental  results 
better  than  Eq.  (24). 
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Much  of  the  experimental  data  is  available  in  terms  of  Hfl,  The 
corresponding  value  of  Rx  can  be  found  from  Eq.  (36)  in  the  same 
manner  in  which  Eq.  (29)  was  obtained  from  Eq.  (27), 


0 


I’p  O.OH  do};  It,'  -  2.3686) 
o  '1  '  hog  Ti^ ' 


After  integration 


Cp  _  ^  ^  0.044 

^  Rx  T'  (log  R/-  1..S)' 


(38) 


(39) 


Therefore 


R 


d  = 


0.044  ' 

dog  '  -  l.S)" 


(40) 


From  this  equation,  the  determination  of  Ri')  is  straightforward  when  Rx' 
is  given;  however,  when  Ri9  is  given,  Rx'  must  be  found  by  some 
iterative  method  such  as  Newton's.  In  a  similar  manner. 


= 


0.044  Uv. 

_ ^ 


(log  It  .  - 
^  1 


(41) 


Once  Rx'  and  Rx.  are  found  to  satisfy  Eqs.  (40)  and  (41),  Cf  and  Cf.  can 
be  found  from  Eqs.  (36)  and  (35)  for  Rf)  =  R)i.  and  a  new  ratio  of  (d  T.f. 

can  be  determined.  Such  a  ratio  is  shown  in  Fig.  3  along  with 
experimental  data  from  Refs.  10,  11,  18,  19,  and  20.  Again,  good 
correlation  is  shown,  indicating  that  Eq.  (36)  is  satisfactory  for  the 
determination  of  Cf. 


DETERMINATION  OF  H, 

As  shown  previously  by  Eq.  (13),  the  transformed  form  factor  can 
be  expressed  in  terms  of  the  incompressible,  adiabatic  form  factor  and 
temperature  ratios.  It  is  further  shown  in  Appendix  B  that  Hi  is  related 
to  Cf.  by 

Hi  =  - (42) 

and  is,  therefore,  a  function  of  the  incompressible  Reynolds  number. 
Thus,  there  remains  the  problem  of  determining  the  relation  between 
the  compressible  and  incompressible  Reynolds  number. 

In  the  transformation  of  Eq.  (3)  into  Eq.  (8)  by  means  of  Eq.  (6) 
and  (7),  no  transformation  of  .v  was  involved.  However,  Ecj.  (8)  can  be 
written  as 
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dX 


+ 


Ue 


dU. 

dX 


(2  +  Hir) 


+ 


dr 

dX 


dx 

dX 


sec  CO 


y+i 


(43) 


where  Ug  is  the  transformed  velocity  (Appendix  A)  and  the  transformation 
dx/dX  must  be  defined.  For  incompressible  flow,  the  momentum 
equation  can  be  written  as 


dxj 


+ 


(2  +  Hi) 


+ 


sec  CO 


(44) 


with  fluid  properties  and  evaluated  at  T^.  Comparison  of  these 
equations  indicates  that  the  transformed  flow  may  be  considered  in¬ 
compressible  if 


dx  (45) 

where  Cf  is  obtained  as  a  function  of  Reynolds  number  based  on  x  and 
Cf.  is  obtained  as  a  function  of  Rx  =  po  Ue  X/po. 

The  original  Stewartson's  transformation  used  the  relation 


sy  -  1 

dX  /  T,  \  2(y-l) 
dx  “  \  To/ 


Culick  and  Hill  (Ref.  21)  using  the  stagnation  temperature  as  a 
reference  temperature  obtained  the  relation 


(46) 


y+  I  N  -  1 

dx  I  tJ 


(47) 


Mager  (Ref.  22)  using  a  still  different  approach  obtained  the  relation 


y+i 

dx  I  To/ 


(48) 


which  reduces  to  Stewartson's  transformation  if  the  viscosity  is  assumed 
proportional  to  the  temperature. 


If  Eqs.  (35)  and  (36)  are  substituted  into  Eq.  (45),  the  result  is, 
since  Rx  =  Hx;  , 


Po  (log  Rx-2.3686)  _  j 
Po  t'e  (logRx-1.5)‘ 


y  +  1 

/  T,  \  2(y-l)  (log  ll„  2.3686) 

V  To/  T'  (log  R  ,  '  -  1.5)' 


P' 

P  Ue 


dRx^  (49) 
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which,  after  integrating,  yields 


Since 


_ Rx _ 

Po  Uedog  Rx  -  1-5)^ 


y+1 

I  To  j 


Tg _ /i "  R  X  ' _ 

T'  p'uodog  Hx'-  l.s)" 


(50) 


1 


11 

Te 


Eq.  (50)  reduces  to 


Rx _ ^  P _ _ 

(log  -  1.5)'  Ho  (log“R^'  -  1.5)' 


(51) 


which  relates  the  incompressible  Reynolds  number  to  the  reference 
Reynolds  number.  Eq.  (51)  cannot  be  solved  directly  for  Itx,  but  the 
use  of  Newton's  approximation  with  Rx  =  ip7po)^^x'  as  the  first 
approximation  yields  as  a  second  approximation 


Rx  = 


p'  ny' 

_ Ifo _ 

l^log  -  2.3686) 


(log  R, 


1.5)' 


-0.8686 


(52) 


which  is  sufficiently  accurate  to  determine  Cfj  from  Eq.  (35).  This 
value  of  Cf.  is  then  used  in  Eq.  (42)  to  determine  Hi,  after  which  FFtr  can 
be  found  from  Eq.  (13)  for  use  in  integrating  Eq.  (10). 


APPLICATION  TO  AXISYMMETRIC  TUNNELS 


In  the  application  of  Eq.  (10)  to  calculate  the  boundary  layer  growth 
of  an  axisymmetric  wind  tunnel,  the  calculations  are  usually  begun  at 
the  throat  where  the  Mach  number  is  unity.  Equation  (10)  can  thus  be 
written 

y+  1 

T^\  2(y-  n 
/ 


y+ 1 

2(y-l)  2+H, 

M 


=  6* 


+ 


y+i 
2  (y  -  1 ) 


M 


2+H„ 


Cf 

2 


sec  oj 


dx 


(53) 
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In  many  cases,  the  momentum  thickness  at  the  throat  d*  can  be  assumed 
to  be  equal  to  zero  without  appreciably  affecting  the  values  calculated 
near  the  end  of  the  nozzle. 


The  Reynolds  number  at  the  throat  must  be  based  upon  a  value  of  x 
which  is  other  than  zero  in  order  to  start  with  a  finite  value  of  Cf.  A 
useful  approximation  for  this  initial  value  of  x*  can  be  obtained  in  the 
manner  suggested  by  Sibulkin  in  Ref.  23.  The  velocity  gradient  at  the 
throat  is  obtained  in  terms  of  the  radius  of  curvature  at  the  throat: 


M  =  1 


y  +  1 


(54) 


It  is  then  assumed  that  the  velocity  gradient  is  constant  from  zero  to 
the  throat  so  that 


Therefore, 


\i  =  I 


X 


* 


y  +  1 
2 


r*  R* 


(55) 


(56) 


The  Reynolds  number  at  each  point  along  the  nozzle  is  thus  based  upon 
a  value  of  x  equal  to  x*  plus  the  distance  from  the  throat  to  the  point. 

In  general,  the  distance  s  along  the  nozzle  contour  maybe  assumed  to 
be  equal  to  the  distance  x  along  the  axis  for  the  purpose  of  determining 
the  Reynolds  number. 


The  above  approximation  can  also  be  used  to  estimate  the  momentum 
thickness  at  the  throat  with  the  additional  assumption  that  the  values 
of  the  reference  temperature  and  the  reference  Reynolds  number  are 
constant  and  equal  to  the  values  at  the  throat.  Equation  (10)  may  then 
be  written 
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and 


du- 


dM 


(59) 


thus,  after  substitution  and  rearranging, 

5  y-3 


(9' 


0.044  (log  R^. 2.3686)  T„  V^*  ^ 


(log  R,.  1.5)’  T*' 


I 


1  .5  +  Htr 
M  dM 


(i  +  V  "’) 


ii  y-  9 
4  (y  -  1 ) 


(60) 


or,  when  y  =  1.4, 

„  0.0694(logR,. 2.3686)  To  R* 

(log  R,^, 1.5)’  T*  J  (1  +  0.2M’)*  (61) 

The  integration  indicated  in  Eq.  (61)  can  be  performed  analytically  for 
values  of  Ibr  which  are  odd  multiples  of  0.5.  For  intermediate  values, 
the  integration  must  be  performed  numerically.  The  results  of  such 
integrations  are  shown  in  Fig,  4, 

For  conical  nozzles  in  which  radial  flow  can  be  assumed,  the 
Mach  number  is  constant  over  each  spherical  segment  of  area  27r  r’/d  +  cos  oj), 
Since  the  throat  area  is  nr*'. 


A  ^  2r’ _ 

A*  (1  +  cos  oj)  r*’ 


M 


1  + 


y-1 

2 

y+1 

2 


M’ 


y+j_ 
2  (y  -  1 ) 


(62) 


or 


r 


+  cos  (U 
2M 


y+1 

2 


M’ 


y+1 

4(y-  1) 


(63) 


from  which 


dr 

dM 


r* 


~r~ 


V 


1  +  cos 
2M’ 


(M= 


-  l)(l  + 


y-1 

2 


5-3y 
ij  4(y-i) 


y+1 

y+  1  4(y-  1) 

2  / 


(64) 


19 


AEDC-TR-59-3 


Then,  since  dr/dx  =  tan  tu 

5  -  3  y 

- hr, - 

y+ 1  ^  4(y-  1)  (65) 

Substituting  Eqs.  (63)  and  (65)  into  Eq.  (10)  yields 


dx  = 


4 


1  -t-  cos  CJ 
2  M" 


(M’  -  l)(l 


1  + 


y+l  n 

“ 

y  +  l  n 

j2(y-l)  ^2  +  H., 

b 

(  Te  > 
1  To  J 

2(y-l)  2  +  H„ 

1  M 

+ 


r* 


(1  +  cos  (u)  CSC  (i) 

1  'j  2(y-i) 


b 

(M^-1)  Cf/2 


1  + 


y-1 

2 


M* 


dM 


(66) 


As  before,  if  the  reference  temperature  and  reference  Reynolds  number 
can  be  assumed  to  be  constant. 


y+i 

/  T,,  x2(y-l) 

I  tJ 


M 


2  +  H„ 


(4t) 


y+l 


1 


2(y-l)  ,,  2  +  H,r 


M 


J 


+ 


r*  (1  +COS  6j)  CSC  <u  Tg  C 

y-*- 1 


dM 


(67) 


and,  if  y  =1.4, 


e 


M 


2  +  H, 


M 


2  +  Hi, 


^  0-00637  (log  R,--  2.3686)  Tor*  cot  ^  rV-i)  M»tr  dM  (68) 

(log  R,  '  -  1.5)*  T'  I  (1  +  0.2  M’)* 

The  integration  indicated  in  Eq.  (68)  can  be  performed  analytically 
when  Htr  is  an  integer  and  numerically  when  Htr  is  not  an  integer.  The 
results  of  such  integrations  are  shown  in  Fig.  5  where  the  integration 
in  each  case  is  performed,  for  convenience,  over  the  interval  from 
1  to  M  .  Obviously, 


/‘’f(M)dM  =  /‘’f(M)dM  -  /‘■f(M)dM 

a  1  1 
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For  the  general  case  of  conical  and  contoured  nozzles,  the  reference 
temperature  and  Reynolds  number  are  not  constant,  and  for  the  con¬ 
toured  nozzle  the  radius  is  not  an  analytic  function  of  the  Mach  number. 
Equation  (53)  must  therefore  be  used  downstream  of  the  throat  while 
Eq.  (60)  or  (61)  can  be  used  to  estimate  the  value  at  the  throat.  In 
many  cases,  the  variation  of  Rx  over  the  length  of  the  nozzle  is  small 
enough  that  a  constant  average  value  of  Hi  can  be  used  and  the  computa¬ 
tions  involved  can  thereby  be  simplified. 

Values  of  boundary-layer  displacement  thickness  calculated  by  the 
method  described  are  compared  with  experimental  values  in  Figs.  6 
and  7  for  Mach  7  and  8  conical  nozzles  and  in  Fig.  8  for  a  Mach  8 
contoured  nozzle.  The  experimental  data  were  obtained  in  the  50-inch- 
diameter  nozzles  of  the  Gas  Dynamics  Facility.  The  agreement  shown 
is  considered  to  be  extremely  good. 


CONCLUDING  REMARKS 


A  method  has  been  developed  for  calculating  the  growth  of  a 
turbulent  boundary  layer  at  hypersonic  Mach  numbers.  Excellent 
agreement  with  experimental  results  from  axisymmetric  nozzles  has 
been  obtained  through  the  application  of  this  method.  The  basis  for  the 
calculations  is  Eq.  (10)  wherein  the  compressible  skin-friction 
coefficient  is  obtained  from  Eq.  (36),  and  the  transformed  form  factor 
is  given  by  Eq.  (13)  in  which  the  incompressible  form  factor  is  related 
by  Eq.  (42)  to  the  incompressible  skin-friction  coefficient  evaluated  at 
the  incompressible  Reynolds  number  given  by  Eq.  (52). 
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APPENDIX  A 

STEWARTSON’S  TRANSFORMATION 


In  the  Stewartson's  transformation,  as  given  in  Ref.  9,  the  distance 
normal  to  the  surface  is  transformed  by  the  relation 


dY  = 


(A-1) 


and  the  velocity  parallel  to  the  surface  is  transformed  by  the  relation 


External  of  the  boundary  layer 


so  that 


(A- 2) 


(A-3) 


(A-4) 


The  definition  of  boundary-layer  momentum  thickness  is 


(A-5) 


where  A  is  the  value  of  y  where  both  velocity  and  temperature  reach 
their  free-stream  values.  Substitution  of  Eqs.  (A-1)  and  (A-4)  into 
Eq.  (A-5)  yields 


6  = 


Po  ap  r^_LL  (i 
Pe  a*  Jo  Ue  ' 


(A-6) 


The  transformed  momentum  thickness  is  defined  as 


0lr  = 


(A-7) 


so  that 


e 


P  0^0 
Pe 


0tr 


(A-8) 
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or 


0  .  (l  *  M.) 


(A-9) 


It  may  be  noted  that  Eq.  (A-7)  has  the  form  for  incompressible  flow, 
although  the  shape  of  the  velocity  profile  is  distorted  by  the  trans¬ 
formation  (see  Ref.  24). 


In  a  similar  manner,  the  boundary-layer  displacement  thickness  is 
defined  by  the  equation 


8* 


(A-10) 


Since  the  pressure  is  assumed  to  be  constant  through  the  boundary 
layer,  Eq.  (A-10)  may  be  written  as 


5* 


(A-ll) 


The  static  temperature  distribution  through  the  boundary  layer  may  be 
expressed  by 


T 


Te 


(> 


(A- 12) 


where  Ts  is  the  local  stagnation  temperature  corresponding  to  the  local 
static  temperature  T.  Substituting  Eqs.  (A-1),  (A-4),  and  ( A-1  2)  into 
Eq.  (A-ll)  yields 


where,  by  definition 


dY 


(A-13) 


(A-14) 


(A-15) 
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When  Ts  is  constant  and  equal  to  To,  Eq.  (A-15)  has  the  form  for  in¬ 
compressible  flow. 

Division  of  Eq.  (A-14)  by  Eq.  (A-8)  yields 


or 


S* 

e 


2 


+ 


y-  I 
2 


(A-16) 


II  =  Mtr  (l  ^  M’)  + 


(A-17) 


Equations  (A-9)  and  (A-17)  were  used  to  transform  Eq.  (3).  Equation 
(A-15)  may  be  written  as 


so  that 


where 


(A-18) 


(A-19) 


(A-20) 


In  order  to  evaluate  Eq.  (A-19),  use  is  made  of  Crocco's  quadratic 
temperature  distribution,  written  as  in  Ref.  8, 


T  _  ^  /  Tyy  —  Tfly,  \  /  U  \  ^  /  T  a  ty _ 1  e  \  /  u 

Ta  Te  \  Ta  /  I  “e  ^  \  Te  ■'  V  '‘e 


{A-21) 


which  assumes  that  the  thermal  boundary  layer  has  the  same  thickness 
as  the  velocity  boundary  layer.  This  would  occur  only  if  the  Prandtl 
number  is  unity.  However,  Eq.  (A-21)  partially  accounts  for  Prandtl 
numbers  other  than  unity  through  the  use  of  the  adiabatic  wall  tempera¬ 
ture  instead  of  the  stagnation  temperature.  Substitution  of  Eq.  (A-21) 
into  Eq.  (A-11),  along  with  Eqs.  (A-1)  and  (A-4)  yields 


(A-22) 
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or 


5* 


(A- 23) 


Division  of  Eq.  (A- 23)  by  (A-8)  yields 


Equations  (A- 24)  and  (A-17)  maybe  combined  so  that 


Hir 


Hi 


1 


(A- 24) 


(A- 25) 


The  variation  of  11  according  to  Eq.  (A-24)  as  a  function  of  Mach 
number  is  compared  in  Fig,  9  with  experimental  results  from  Refs.  10 
and  25  and  in  Fig.  10  with  the  tabulated  results  from  Ref.  26.  The 
correlation  is  considered  to  be  very  good  and  justifies  the  derivation 
described. 
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APPENDIX  B 

INCOMPRESSIBLE  SKIN-FRICTION  COEFFICIENTS 


A  great  many  empirical  equations  have  been  developed  for  the 
variation  of  the  local  and  mean  skin-friction  coefficients  with  Reynolds 
number  for  incompressible  flow.  Probably  the  most  simple  to  use  are 
the  power-law  equations 


2 


(B-1) 


(B-2) 


where 


CPi  _  N  w  Cf. 

2  K,.  ~  N  2 

^  1 


(B-3) 


/3 


(B-4) 


For  a  limited  range  of  Reynolds  numbers,  values  of  a,  /3,  and  .'N  can  be 
considered  to  be  constant;  but,  for  a  large  range  of  Reynolds  numbers, 
their  values  must  also  be  functions  of  the  Reynolds  number. 


Locke  (Ref.  27)  investigated  a  number  of  the  empirical  equations 
and  found  that  Schoenherr's  equation 


0.242 

V^F  i 


log  (CF;  llxj) 


(B-5) 


gave  a  good  correlation  with  a  very  large  amount  of  experimental  data 
over  the  wide  range  of  Reynolds  numbers  from  10'  to  10”.  Other  forms 
of  this  equation  are 


0.242 


log  (2  R0.) 


(B-6) 
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and 


‘i  log  (  2  R(9.  )  [  log  (  2  R|3,  )  +  0.8686  ]  lB-7) 

Although  Schoenherr's  equations  correlate  well  with  experimental 
results,  they  are  difficult  to  use  because  neither  Cf.  or  Cfj  can  be  found 

directly  as  a  function  of  Rx;  .  Tables  can  be  constructed  for  use  when  the 

computing  is  done  manually;  however,  more  direct  equations  are  desired 
for  use  with  automatic  computers.  Any  set  of  equations  must  satisfy  the 
relationship 


or 


Cf.  =  Cp.  +  R, 


d  Cl 


d  R, 


(B-8) 


(B-9) 


After  an  investigation  of  several  types  of  equations,  the  equations 


and 


0.088 

(log  R^.  -  1.5)" 


(B-10) 


Cf. 


0.088  (log  R^.  -  2.3686  ) 

(log  R,,  -  1.5)’ 

1 


(B-11) 


were  found  to  correlate  with  experimental  data  as  well  as  Schoenherr's 
equations  over  the  range  of  Reynolds  numbers  from  10'  to  10’  .  A 
comparison  of  values  of  Cp;  calculated  from  Eq.  (B-10)  with  those  from 

Eq,  (B-5)  is  shown  in  Fig.  11  together  with  the  "ideal"  values  from 
Ref.  28  and  29.  A  similar  comparison  of  values  of  Cf.  is  shown  in 

Fig.  12  together  with  experimental  values  of  Dhawan  (Ref.  30),  Schultz- 
Grunow  (Ref.  31),  and  Kempf  (Ref.  3  2).  These  data  are  recognized  as 
being  among  the  most  accurately  determined  values  available.  A 
further  comparison  of  R(9.  as  a  function  of  R*;  is  made  in  Fig.  13  with 

the  experimental  values  of  Wieghardt  (Ref.  33).  All  of  these  com¬ 
parisons  indicate  the  validity  of  Eqs.  (B-10)  and  (B-11). 


A  critical  examination  of  the  curve  of  Eq.  (B-10)  shows  that  it 
crosses  Schoenherr's  curve  at  a  Reynolds  number  of  about  10'  and  again 
at  about  5  x  10’  .  Moreover,  it  goes  to  infinity  when  log  Rx;  =1.5  or 

Rxj  =32.  Furthermore,  the  curve  of  Eq.  (B-11)  has  zero  slope  when 
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log  Rxj  =  2.  3686  or  Rx;  =  234.  Both  of  these  conditions  are  well  below 

the  range  for  turbulent  boundary  layers  and  are,  therefore,  of 
academic  interest  only.  The  values  of  the  constants,  0.  088  and  1.  5, 
could  have  been  chosen  to  agree  better  with  Coles'  values  but  the 
experimental  values  do  not  warrant  such  a  selection. 


The  incompressible  form  factor  Hj  is  also  a  function  of  Reynolds 
number  and  is  shown  in  Ref.  29  to  be  related  to  the  skin-friction 
coefficient  in  the  manner. 


Hi 


1 _ 

1  -  7  V  Ci./2 


(B-12) 


where  the  constant  7  (although  different  from  that  of  Ref.  29)  is  chosen 
for  correlation  with  experimental  data.  Such  correlations  are  shown 
in  Figs.  14  and  15.  Again  the  correlations  are  quite  satisfactory. 

For  comparative  purposes,  values  of  Cp.,  Cfj,  R^.  ,  and  Hi  are 

listed  in  Table  1  for  various  values  of  Rx.  from  10*  to  10’ .  Also 

listed  are  the  values  of  N  obtained  from  Eq.  (B-3)  which  show  that  the 
use  of  Eqs.  (B-1)  and  (B-2)  must  be  limited  to  small  ranges  of 
Reynolds  number. 
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TABLE  1 

INCOMPRESSIBLE  SKIN  FRICTION  VALUES 


10^.  .0051j01  .00718U  359.2  1.5726  3.03 

2  X  lo5  .00ii699  .006091  609.1  1.5135  3.38 

3  X  10^  .Ooa3U8  .00556U  831.5  I.li81i5  3.58 

h  X  lo5  .OOU122  .005230  1016.0  1.U658  3.72 

5  X  105  .003959  .001991  1218.0  l.ii523  3.83 

-  6  X  105  .003832  .00h808  lliii2.0  l.l')il8  3.93 

7  X  lOj  .003731  .OOL662  1632.0  l.U33ii  li.OO 

8  X  105  .OO36UU  .00^539  1816.0  I.I126I  h.07 

9  X  105  .003571  .OOIU35  1996.0  l.ij200  ii.13 

10®,  .003507  .001^316  2173.0  l.iill47  1.18 

2  X  10®  .003127  .003818  3818.0  1.3827  U.53 

3  X  106  .002986  '  .003596  53^'.0  1.3707  h.73 

h  X  lof  .002805  .003381  67^.0  1.3553  U.87 

5  X  10®  .002712  .003256  8139.0  1.3U73  1^.99 

6  X  10°  .002639  .003159  9U76.O  1.3U0  5.08 

7  X  10®  .002580  .003080  10780.0  1.3359  5.15 

8  X  10®  .002530  .003011;  12060.0  1.3315  5.22 

9  X  106  .OO2U87  .002958  13310.0  1.3277  i  5.28 

107  .002l;50  .002909  ll;550.0  1.32l;5  5.33 

2  X  107  .002223  .002615  26150.0  l,30hh  5.68 

3  X  lo7  .002105  .002163  36950.0  1.2938  5.88 

h  X  107  .  002027  .  002363  1j7270.0  1.2868  6.03 

5  X  107  .001969  .002290  57250.0  1.2815  6.II; 

6  X  107  .OOI92U  .002233  66980.0  1.2773  6.23 

7  X  107  .001887  .002186  76500.0  1.2739  6.31 

8  X  107  .001855  .00211;6  85850.0  1.2709  6.37 

9  X  1Q7  .001828  .002112  95060.0  1.2681;  6.1;3 

10°  .001805  .002083  lOhlOO.O  1.2663  6.1;8 

2  X  108  .001660  .001903  190300.0  1.2526  6.83 

1;  X  10®  .001531  .001715  31;8900.0  1.2l;02  7.18 

6  X  108  .001i;63  .001661  1;981;00.0  1.2335  7.38 

8  X  108  .0011;17  .001606  61;2300.0  1.2290  7.52 

109  .001383  .001561;  782200.0  1.2256  7.63 
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O  MEASURED  DISPLACEMENT  THICKNESS 
-  PRESENT  THEORY 


M  =  7.07 
Tw=  520®  R 
To=  II20®R 
po=  390  PSIA 
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Fig.  6  Comparison  of  Calculated  Boundary-Layer  Thickness  with  Values  Measured  in  a  Mach  7 
50-inch-Diameter  Conical  Nozzle 
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Fig.  7 


Comparison  of  Calculated  Boundary-Layer  Thickness  with  Values  Measured  in  a  Mach  8, 
50-inch-Diameter  Conical  Nozzle 
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BOUNDARY-LAYER  FORM  FACTOR 
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Fig.  9  Correlation  of  Calculated  Values  of  H  with  Experimental  Values 
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D  REF.  34  SERIES  T 
0  REF.  3  4  SERIES  E 
0  REF  3 3, 33.0  METERS/SEC 
0  REF.  3  3, 1 7. 8 METERS/SEC 
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